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ABSTRACT 



O 
U_ 

r^ I Aims. Within the framework of the HERM33ES Key Project, using the high resolution and sensitivity of the Herschel 

O . photometric data, we study the compact emission in the Local Group spiral galaxy M 33 to investigate the nature of 

I ■ the compact SPIRE emission sources. We extracted a catalogue of sources at 250 /im in order to investigate the nature 

O ' of this compact emission. Taking advantage of the unprecedented Herschel resolution at these wavelengths, we also 

focus on a more precise study of some striking Ha shells in the northern part of the galaxy. 
C/) , Methods. We present a catalogue of 159 compact emission sources in M 33 identified by SExtractor in the 250 /xm SPIRE 

Cu . band that is the one that provides the best spatial resolution. We also measured fluxes at 24 /xm and Ha for those 159 

extracted sources. The morphological study of the shells also benefits from a mult i wavelength approach including Ha, 
far- ultraviolet from GALEX, and infrared from both Spitzer IRAC 8 /im and MIPS 24 fim in order to make comparisons. 
Results. For the 159 compact sources selected at 250 /xm, we find a very strong Pearson correlation coefficient with the 
MIPS 24 //m emission (r24 = 0.94) and a rather strong correlation with the Ha emission, although with more scatter 
{vHa = 0.83). The morphological study of the Ha shells shows a displacement between far- ultraviolet. Ha, and the 
SPIRE bands. The cool dust emission from SPIRE clearly delineates the Ha shell structures. 

Conclusions. The very strong link between the 250 /xm compact emission and the 24 /xm and Ha emissions, by recovering 
the star formation rate from standard recipes for H ii regions, allows us to provide star formation rate calibrations based 
l/^ ■ on the 250 /xm compact emission alone. The different locations of the Ha and far- ultraviolet emissions with respect to 

^^ ' the SPIRE cool dust emission leads to a dynamical age of a few Myr for the Ha shells and the associated cool dust. 
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1. Introduction mation (SF). Indeed, the distance, 840 kpc (Freedman et al. 

;J] ^^^. , . , p 1 r 1 . ^ . 1991), and inclination, 56 degrees (Regan & Vogel 1994), of 

=^.!^'.* .*,.*; 0^"'''^. ° ^ "". °P'''' ^''^f ""^ P'°Jf* M 33 allow us to reach an unprecedented spatial resolution 

■ - _ HerschelM 33 extended survey (HERM33ES) , we are study- j^ ^j^^ far-infrared of an external spiral galaxy and investi- 

mg the galaxy M 33 to understand the origin of various di- ^^^ ^j^^ phenomena in the interstellar medium which lead 

agnostic hues in heating and coohng and other processes ^^ ^j^^ gp ^-^^ Herschel (Pilbratt et al. 2010), we are able 

m the interstellar medium (ISM . We refer the reader to ^^ ^^^^^^^ 35 ^09 ^^^ ^70 structures of the cool 

Kramer et al. (2010) for more details about the overall ^^^^ -^ ^^^ gPIRE 250, 350, and 500 ^^m bands, respec- 

^^^^^f^l P'°f ^* ?°^J;'' ''' '^''^'^1 ^°' "" ^''* presentation of ^-^^^^^ ^j^j^j^ -^ ^^^ ^ypj^^l ^p^^-^j ^^^j^ ^f gj^^^ jj „ ^^gj^^^ 

the PACS and SPIRE maps of the entire galaxy, together ^^^ ^j^^-^ complexes. In this article, we follow a two-fold ap- 

with spatially averaged spectral energy distributions. p^^^^j^ -^ ^^^^^ ^^ ^^^^^ ^^^ ^^^^^^ ^j ^1^^ ^^^p^^^ 250 ^im 

The local group late-type spiral galaxy -M 33- is an gj^ission and take advantage of the resolution and sensitiv- 

ideal target for studying the detailed processes of star for- -^y ^^^^j^^j ^^ Herschel to resolve the cool dust emission 



* Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consor- 
tia and with important participation from NASA. 



associated with H ll regions. 
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Fig. 1. SPIRE maps of M33 at 250 /im (logarithmic scale). 
The catalogue of extracted sources is superimposed. The 
ellipses are indicative radii and are meant to show the po- 
sition of the sources, they are more extended than where 
the photometry is actually performed. The white rectangle 
delimits the zone detailed in Fig. 3. 



2. Data 

2.1. Herschel photometric data 

During the photometric observations of the HERM33ES key 
project, we mapped M33 with PACS and SPIRE instru- 
ments in parallel mode in two orthogonal directions, inte- 
grating 6.3 hours (two AORs). In this letter, we particularly 
focus on the nature of the compact emission in the SPIRE 
bands. The reduction includes all basic pipeline steps plus 
a robust polynomial fit for the baseline, which improves 
the images by reducing scan line structures (more details 
about the SPIRE observations are given in Kramer et al. 
2010). The r.m.s. noise levels of the SPIRE maps are 14.1, 

9.2, and 8mJy/beam, at 250, 350, 500 /im. The mean mea- 
sured spatial resolutions for Gaussian beams are 19'.'7±1'.'9, 
27:^0 ± 4'.'1, 37V2 ± 2V3 for 250, 350, and 500 /im, respec- 
tively. The overall absolute calibration accuracy is esti- 
mated as within ±15%, following the SPIRE instrument 
paper (Griffin et al. 2010). 



2.2. Ancillary multiwavelength data 

To perform comparisons with the Herschel far-infrared 
emissions, we use the following mult i- wavelength ancillary 
data: Ha narrow-line image (Hoopes & Walterbos 2000) 
to trace the current unobscured SF (<10Myr), GALEX 
far- ultraviolet (Gil de Paz et al. 2007) to trace the recent 
unobscured SF (<100Myr), and Spitzer IRAC 8/im and 
MIPS 24 /im (Verley et al. 2007, 2009, 2010) infrared data 



to trace the obscured (probably recent) SF. In the Ra 
map, the [Nil] contamination is less than 5%, and we have 
corrected the flux for the foreground Galactic extinction, 
Aua =0.17 (Relaho & Kennicutt 2009). We consider typi- 
cal uncertainties to be better than 15% for the MIPS 24 /im 
(Verley et al. 2009) and the Ha photometric measurements 
(Hoopes & Walterbos 2000). 

3. Nature of the compact SPIRE 250 /im emission 

In this section, we investigate the nature of the compact 
SPIRE 250 /im emission, in a statistical manner. To create a 
catalogue of compact emission sources in the SPIRE 250 /im 
band, we use the SExtractor software (Bertin & Arnouts 
). The photometry of the sources is computed using 
the parameter FLUX_IS0 given by SExtractor, which uses 
isophotal photometry (sum of all the pixels above a thresh- 
old given by the lowest isophot: 16 times the background 
r.m.s.). Of course, the number, position, size, and photom- 
etry of the obtained sources can vary widely by choosing 
other input parameters in order to pursue a given scientific 
goal, but the above parameters are well-suited to recover 
as much of the compact 250 /im emission regions as we can, 
and lead to a final catalogue of 159 sources. The catalogue 
of extracted sources is overplotted in Fig. 1. 

Interestingly, we detect almost no sources in the inter- 
arm regions. Since most of the extracted objects are along 
the spiral pattern of the galaxy, we believe many of the 
sources may be directly linked to SF. This suggests that 
the SPIRE 250 /im compact emission could be a reliable 
SF tracer in the vicinity of H ll regions. Special care must 
be taken for the most diffuse part of the 250 /im emission, 
which may be powered by a more general interstellar radia- 
tion field and may not be directly linked to current or recent 
SF (see for instance Israel et al. 1996). Detailed studies of 
the reliability of the diffuse PACS and SPIRE emission as 
SF tracers will be presented and discussed in related pa- 
pers (Boquien et al. 2010, Boquien et al., Calzetti et al., in 
prep.). For the present letter, we concentrate our prelimi- 
nary work on the 250 /im compact emission and compare its 
properties with standard SF tracers such as the Ha emis- 
sion line and the 24 /im compact emission linked to Hll 
regions (Calzetti et al. 2005, 2007). 

To do so, we performed photometry on the Ha and the 
Spitzer MIPS 24 /im images, using SExtractor to reproduce 
photometry on different images with the same centres and 
apertures. Thus, after having degraded (using a Gaussian 
function) the images to the 250 /im SPIRE resolution, we 
obtained the Ha and 24 /im fluxes with the same input 
parameters as the ones used in the SPIRE 250 /im im- 
age. In the upper panel of Fig. 2, we show the Ha and 
24 /im luminosities as a function of the SPIRE 250 /im lu- 
minosity for the final 159 sources. Remarkable correlations 
appear between the Ha and 24 /im luminosities and the 
250 /im luminosity of these sources, confirming that the 
cool dust compact emission, as traced by the bright SPIRE 
250 /im emission, is closely linked to SF. The 24 /im emis- 
sion presents a very high Pearson correlation coefficient 
(r24 = 0.94, using the logarithmic values of the luminosi- 
ties, as displayed in Fig. 2), while, although rather strongly 
correlated {rua = 0.83) with the 250 /im emission, the Ha 
presents more scatter. That the 250 /im emission correlates 
better with the 24 /im emission is not surprising because 
both are linked to dust, either cool for 250 /im or warm for 
the 24 /im, while the Ha emission directly reflects the ion- 
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ising photons of hot OB stars. Part of the scatter of the Ha 
observed luminosity can also be due to emission extincted 
by dust. A certain scatter between L(250/im) and SFR is 
therefore expected because of variations in dust tempera- 
ture, metallicity, maybe a varying initial mass function, and 
although minimised in our analysis, possibly heating due to 
the interstellar radiation field. 

As the Ha and 24 /im compact emissions are standard 
star formation rate (SFR) tracers, we can use the calibra- 
tions given by Calzetti et al. (2007, 2010) for Hll regions 
to obtain the SFR of these sources. To calibrate the SFR 
from the 24 /im emission, we use Eq. 13 in Calzetti et al. 
(2010) while for the combined Ha+24/im SFR calibration, 
we use their Eq. 16. We can therefore obtain the calibra- 
tion to recover the SFR from the 250 /im emission alone, by 
matching the existing Ha+24/im and 24 /im SFR for Hll 
regions. To recover the Ha+24 /im SFR, our best fit leads 
to 



SFR [Moyr-^] = 8.71 x 10-''^L(250/im) 



\1.03 



(1) 



where 1/(250 /im) is in erg s ^. To recover the SFR(24), we 
need the following calibration: 

(2) 



SFR [Mo yr-^] = 3.47 x 10-'^'^L(250 /im)^-^^ , 

also with L(250/im) in erg s~^. The uncertainties are 0.04 
and 0.03 for the exponents in Eqs. 1 and 2, respectively, 
while the calibration constants have uncertainties of 4.0 and 
2.7%, respectively. The two equations show rather similar 
exponents, compatible to unity within la, and reflect the 
almost linear behaviour of the compact 250 /im emission 
with respect to SFR. These recipes would then be valid to 
obtain SFR in H ll regions using the compact 250 /im emis- 
sion for 10^^ < L(250/im) < lO^^'^ergs"^, approximately. 

4. SPIRE emission distributions for H ii regions 

The high resolution of the Herschel data for M 33 allows us 
to perform a deeper comparison of the emission distribu- 
tions within the star- forming regions. Following the work 
in Relaho & Kennicutt (2009) and Relaho et al. (2010), we 
focus on SPIRE emission in the interior of the most lu- 
minous Hll regions. Herschel wavelengths are particularly 
well-suited for this study because they trace the cool dust 
components that play a primordial role in H ll regions. From 
inspection of the two most luminous Hll regions in M33, 
NGC 604 and NGC 595, we see that in general the SPIRE- 
bands are displaced from the Ha and FUV emission in the 
interior of the regions and have a more diffuse component 
extended towards the outer parts. Similar displacement be- 
tween the infrared and FUV emissions has already been 
noticed by Calzetti et al. (2005) for about half of the stud- 
ied Hll regions in M51a. As Ha corresponds to the last 
- lOMyr of SF and FUV to the last - 100 Myr, we can 
propose two possible explanations for the displacement: 1) 
as H2 forms on the dust grains, then everything (dust and 
gas) collapses to form stars, and it takes between 10 and 
100 Myr for all the reservoir of dust in a given region to be 
consumed this way; 2) the UV field emitted mainly by OB 
stars and other non-ionising stars can be so strong that it 
pushes away the dust from the star-forming region in less 
than - 100 Myr. 

Support for the second explanation is the emission dis- 
tribution for the most prominent shells in the north part 
of M33 presented in Fig. 3. We compare in this figure the 
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Fig. 2. Upper panel: Ha and 24 /im emissions as a function 
of the 250 /im SPIRE emission for the 159 extracted sources. 
Bottom panel: Star formation rate from 24 /im and from 
Ha+24/im emissions as a function of the 250 /im SPIRE 
emission. Typical uncertainties of 15% are shown in the 
bottom right corner of each panel. 



emissions at Ha, FUV, 24 /im, and SPIRE bands in a set 
of Ha-shells located in the north part of the galaxy. The 
location of these shells within the galaxy (white rectangle 
in Fig. 1) is in a region where no significant 8/im or 24 /im 
emission has been detected, but where diffuse emission is 
observed in the SPIRE bands. The PACS bands (100 /im 
and 160 /im) were also examined, and we found some dif- 
fuse emission in the 160 /im but no significant emission in 
the 100 /im band. As shown in Fig. 3, the shells have strong 
FUV emission in their centres. The northern shell seems to 
have a knot where the emissions at FUV and Ha coincide, 
but the FUV is very nicely located in the centre for the 
bigger middle shell, suggesting that the shell could have 
been created by the stellar winds (SW) or supernova (SN) 
explosions coming from the stars in the cluster emitting 
at FUV. Using the Ha image, the radius of the central 
Ha shell is estimated to be ^150 pc. Assuming an expan- 
sion velocity of '-«^50kms~^ for the Ha shell consistent with 
previous observations of Ha expanding shells in Hll re- 
gions (Relaho & Beckman 2005), we derive a kinematical 
age for the central shell of ^3 Myr. Dynamical ages of a few 
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Fig. 3. Top-left: Continuum-subtracted Ha image of a set of large shells in the outer north part of M33. Contours are 
overlaid in this image to better enhance the shell features (levels are at 20, 49, 122, 300 of Ha emission measure) and 
repeated in the other images for comparison. From top to bottom and from left to right: FUV (res. 4^2), 24 /im (res. 
b'!l), 250 /im (res. 19V7), 350 /im (res. 27V0), and 500 /im (37^2) images. 



Myr are consistent with ages estimated for shells in dwarf 
galaxies (Martin 1998). It is remarkable that no emission at 
24 /im is seen in the central and southern shells, but there 
is emission from cool dust in the three SPIRE bands. We 
suggest that, during the time the shell is forming, the dust 
is mixed with ionised gas inside the H ll region, a dust frac- 
tion is heated by the ionising radiation and emits at 24 /im, 
while the rest is cool and emits in the SPIRE bands. This 
could be the case for the northern shell in the process of 
formation. Later, when the shell is created by the SW and 
SNe, we observe FUV emission in the centre. Ha emission 
from the shell, and cool dust outlining the Ha structure 
(see the central shell where the 250 /im emission clearly de- 
lineates the Ha shell). This would be the first evidence that 
the dust participates in the kinematics of the gas inside the 
H II regions and that the impact of the SW and SNe into the 
ISM also affects the interstellar dust. The behaviour shown 
in Fig. 3 is also seen in a significant number of shells in 
the outer parts of the galaxy, and it would be very interest- 
ing to relate it with the existence of Hi bubbles observed 
in other galaxies (Weisz et al. 2009). A multi- wavelength 
comparison in a larger sample of shells would be needed to 
study this issue further. 

5. Conclusions 

For the Local Group spiral galaxy M 33, using the unprece- 
dented resolution and sensitivity of the Herschel SPIRE 
photometric data, we focus on the compact emission and 
conclude the following. 

— We created a catalogue of 159 sources from the SPIRE 
250 /im emission using a high detection threshold to fo- 
cus on the compact emission of the cool dust; the com- 
pact emission sources associated to the cool dust follow 
the spiral pattern of the galaxy. 

— Very high Pearson correlation coefficients, r24 = 0.94 for 
the 24 /im emission and rHa = 0.83 for the Ha emission, 
confirm that the 159 sources are closely linked to SF. 

— By using standard (Ha+24/im and 24 /im alone) SFR 
recipes for H ll regions provided by Calzetti et al. (2007, 
2010), we are able to calibrate the 250 /im compact emis- 



sion as an SFR tracer and provide conversion factors 
that can be used in further studies (Eqs. 1 and 2). 

— The morphological study of a set of three Ha shells 
shows that there is a displacement between far- 
ultraviolet and the SPIRE bands, while the Ha struc- 
ture is in general much more coincident with the cool 
dust. 

— The different locations of the Ha and far-ultraviolet 
emissions with respect to the SPIRE cool dust emis- 
sions leads to a dynamical age of a few Myr for a set of 
Ha shells and the associated cool dust. 
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